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We report a measurement of the inclusive semileptonic branching fraction of the B s meson using 
data collected with the BABAR detector in the center-of-mass (CM) energy region above the T(4S) 
resonance. We use the inclusive yield of <f> mesons and the (f> yield in association with a high- 
momentum lepton to perform a simultaneous measurement of the semileptonic branching fraction 
and the production rate of B s mesons relative to all B mesons as a function of CM energy. The 
inclusive semileptonic branching fraction of the B 3 meson is determined to be B(B 3 — > ivX) = 
9-5i^o( s t a t)^i'9( s yst)%, where t indicates the average of e and /i. 



PACS numbers: 14.40. Nd 13.20.He 

Semileptonic decays of heavy-flavored hadrons serve 
as a powerful probe of the electrowcak and strong inter- 
actions and are essential to determinations of Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements (see, for 
example, "Determination of V c b and V u b" in Ref. [1]). 
The inclusive semileptonic branching fractions of the Bd 
and B u mesons are measured to high precision by ex- 
periments operating at the T(4S') resonance, which de- 
cays almost exclusively to BB pairs (here and through- 
out this note, BB refers to B^Bj, and B U B U ). How- 
ever, lacking an analogous production mechanism, in- 
formation on branching fractions of the B s meson re- 
mains scarce nearly two decades after its first observa- 
tion [1]. Here we report a measurement of the inclusive 
semileptonic branching fraction of the B s meson using 
data collected with the BABAR detector at the PEP-II 
asymmetric-energy electron-positron collider, located at 
the SLAC National Accelerator Laboratory. The data 



were collected in a scan of center-of-mass (CM) energies 
above the T(4S') resonance, including the region near the 
B S B S threshold. As <f> mesons are particularly abundant 
in B s decays due to the CKM-favored B s — > D s transi- 
tion, the inclusive production rate of <f> mesons and the 
rate of mesons produced in association with a high 
momentum electron or muon can be used to simultane- 
ously determine the B s semileptonic branching fraction 
and the B s production fraction as a function of the CM 
energy E C m- 

The energy scan data correspond to an integrated lu- 
minosity of 4.25 fb" 1 collected in 2008 in 5MeV steps 
in the range 10.54 GcV < E C m < 11.2 GcV. In a previ- 
ous study [2] , we presented a measurement of the inclu- 
sive b quark production cross section Rb = a(e + e~ — > 
bb)/o-°(e + e~ — > m this energy range, using this 

same data sample (tr° is the zeroth-order QED cross- 
section). In the present study, we also make use of 
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18.55 fb 1 of data collected in 2007 at the peak of the 
T(4S) resonance, and 7.89 ftT 1 collected 40 MeV be- 
low the T(45), to evaluate backgrounds from continuum 
(e + e~ — > qq, q = u,d,s,c quark production) and BB 
events. We choose below-resonance data for which de- 
tector conditions most closely resemble those of the scan, 
and on-resonancc data corresponding to roughly twice 
the luminosity of the below-resonance sample. The sizes 
of these samples are sufficient to reduce the correspond- 
ing systematic uncertainties below those associated with 
irreducible sources. 

The BaBar detector is described in detail elsewhere [3] . 
The tracking system is composed of a five-layer sili- 
con vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) in a 1.5-Tcsla axial magnetic field. The SVT 
provides a precise determination of the track parame- 
ters near the interaction point and standalone tracking 
for charged particle transverse momenta (p t ) down to 
50MeV/c. The DCH provides a 98% efficient measure- 
ment of charged particles with p t > 500MeV/c. The 
p t resolution is <T Pt /p t = (0.13 • p t + 0.45)%. Hadron 
and muon identification in BABAR is achieved by using a 
likelihood-based algorithm exploiting specific ionization 
measured in the SVT and the DCH in combination with 
information from an instrumented magnetic-flux return 
and the Cherenkov angle obtained from the detector of 
internally reflected Cherenkov light. Electron identifica- 
tion is provided by a combination of tracking and in- 
formation from the CsI(Tl) electromagnetic calorimeter, 
which also serves to measure photon energies. For the 
evaluation of event reconstruction efficiencies across the 
scan range, simulated samples of e + e~ — > H + jJi~, con- 
tinuum, and e+e~ — > B q *^ B q *\ q — u,d,s events, cre- 
ated with the kk2f [4], JETSET [5], and EvtGen [6] 
event generators, respectively, are processed through a 
Geant4 [7] simulation of the BaBAR detector. 

For this measurement, we present the scan data as a 
function of Ecm in bins of 15 MeV. In each bin we 
measure the number of BB-\ike events (defined below), 
the number of such events containing a <f> meson, and the 
number of events in which the <fi meson is accompanied by 
a charged lepton candidate. The results are normalized 
to the number of e + e~ — > events in the same en- 

ergy bin so that the luminosity dependence in each bin is 
removed. These three measurements are used to extract 
the fractional number of B S B S events and the semilep- 
tonic branching fraction B(B S — > tvX). The procedure 
is described in detail below. 

To suppress QED background, events are preselected 
with a multihadronic event filter optimized to select BB 
and B S B S events. The filter requires a minimum number 
of charged tracks in the event (3), a minimum total event 
energy (4.5 GeV), a well-identified primary vertex near 
the expected collision point, and a maximum value of 
the ratio of the second to zeroth Fox Wolfram moments 
[8] (i?2 < 0.2) calculated in the CM frame using both 



charged tracks and energy depositions in the calorimeter, 
where the latter are required not to be associated with a 
track. 

A different preselection is used to identify muon pair 
events. Events passing this selection must have at least 
two tracks. The two highest momentum tracks arc re- 
quired to be back-to-back in the CM frame to within 
10 degrees, appear at large angles to the beam axis 
(| cos#cm| < 0.7486), and have an invariant mass greater 
than 7.5 GeV/ c 2 . In addition, we require that less than 
1 GeV be deposited in the electromagnetic calorimeter. 
This selection is 43% efficient for simulated /x 4 ' pT events 
while rejecting virtually all continuum events. 

Candidate <j> mesons are reconstructed in the (f> — > 
K + K~ decay mode, by forming pairs of oppositely 
charged tracks that are consistent with the kaon hy- 
pothesis. In each event, the (f> candidate with the best- 
identified K ± daughters is selected by assigning a weight 
to each based on the particle identification criteria. 
The 4> candidate with the largest sum of kaon weights is 
selected. The invariant mass distribution of these can- 
didates is used to determine the 4> yield in a given Eqm 
bin using a maximum likelihood fit. Events containing cf> 
candidates and an electron or muon candidate with a CM 
momentum exceeding 900MeV/c are used to determine 
the yield of events with both a <fi and a lepton (0-lepton 
events) . The requirement on the lepton momentum sup- 
presses background from semileptonic charm decays. 

Figure 1 shows, as an example, the K + K~ invariant 
mass distribution for (a) all (f> candidates, and (b) <p- 
lcpton candidates, in the energy bin 10.8275 < Ecm < 
10.8425 GeV. These mass distributions are fit to the func- 
tion 

f(M; N, b, c) = NV(m KK ; m , r , a) 

+ Nc(l + bm KK Jl-(^) 2 , 

V \ m KK J 

(1) 

with ttik the world-average mass value [1] of the K ± . 
V(m,KK', m^, Tfi, <j) is a Voigt profile (the convolution of 
a Brcit-Wigner function l/((m KK — m^) 2 + r^, 2 /4) with 
a Gaussian resolution function) normalized to unity, so 
that N is the number of events in the peak. We fix the 
mean (m^) and Breit-Wigner width (r^) to the world av- 
erage values of the <f> mass and natural width [1], and the 
width of the Gaussian resolution (er) by first performing 
all of the <j> fits with the parameter left free, then fixing it 
to the weighted mean of all of the values obtained across 
the scan. The value in data determined by this method 
is a = 1.61 ± 0.04(stat) MeV/c 2 . The combinatoric back- 
ground is modeled as the product of a linear term and a 
threshold cutoff function parameterized by the slope of 
the linear term (b) and a relative scaling (c). 

To determine the <p and (/(-lepton yields from B decays 
in each Ecm bin, the contribution of continuum events is 
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FIG. 1. Invariant mass distribution of — >■ K + K candi- 
dates in the energy bin 10.8275 GeV < E C m < 10.8425 GeV: 
(a) inclusive 4> candidates; (b) 0-lepton candidates. The back- 
ground shape is shown by the dashed curve and the total fit 
by the solid curve. 
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FIG. 2. Relative (a) event, (b) (f>, and (c) <j!>-lepton yields, 
normalized to the yields. Corrections for detector effi- 

ciency have not been applied. The dotted vertical line indi- 
cates the B s production threshold. 



subtracted. This is achieved by using the data collected 
below the T(4S f ) described above. The event, 0, and <p- 
lepton yields are measured in this dataset following the 
same procedures described above. These yields are cor- 
rected for the energy dependence of the reconstruction 
efficiencies and are then subtracted from the scan yields 
in each Eqm bin. This procedure neglects the different 
energy dependence of a small component of the hadronic 
and dimuon cross sections, primarily due to the presence 
of initial state radiative (ISR) e+e" -> tT(15, 2S, 3S) 
and two photon e + e~ — > e + e~7*7* — > e + e~Xh events, 
which do not scale according to 1/E^ M . The effect of 
these contributions is to introduce a small energy depen- 
dence on the amount to be subtracted from each bin. The 
average size of this effect is estimated to be less than 2% 
of the below-resonance event yield. The impact on the 
result is taken as a systematic uncertainty. 

The normalized event, (f>, and 0-lepton yields after the 
continuum subtraction are presented in Fig. 2. These 
three quantities, denoted Ch, and C^i respectively, 
can be expressed in terms of contributions from events 



containing B u d and B s events, the cross section ratio 
Rb = <r{e + e~ — > BqBq)/a IM + IM -, and the related 

q={u,d,s} 

reconstruction efficiencies, as follows: 

C h = R B [f s e s h + (1 - f s )e h ] (2) 
= Rb [f s e%P(B s B s -> (f>X) 

+(1 - f s )e^P{BB -> <t>X)] (3) 
C#= Rb [f s e s ^P(B s B s -> <t>lX) 

+(1 - f s )e H P{BB -> tfAX)] (4) 

(with energy dependence implicit in all terms here and 
elsewhere), where 

f EE NB - (5) 

Js N Bll +N Bd +N Bs { ' 

and exi^'x) ^ s the efficiency for a B u ^ (B s ) pair to con- 
tribute to the event, 4> or </>-lepton yield. The efficien- 
cies are estimated from simulation, while P(BB —> 0X) 
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and P(BB — > <fi£X), which are the probabilities that a <f> 
or a 0-lepton combination is produced in an event with 
a BB pair, are measured using the T(45) data sample 
described above. Specifically, we determine the <fr and <p- 
lepton yields in the the T(45) data. We then apply Eqs. 
(2), (3), and (4) with f s = to extract e^P{BB -> 4>X) 
and t^iPiBB — > <fi£X). Simulations are used to extrapo- 
late the values of the efficiencies to other energies. 



— 0.2 
0.15 
0.1 
0.05 


-0.05 
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FIG. 3. Results for the fraction / s as a function of Ecm- 
The inner error bars show the statistical uncertainties and 
the outer error bars the statistical and systematic uncertain- 
ties added in quadrature. The dotted line denotes the B a 
threshold. 



The remaining unknown quantities of interest are the 
probabilities P(B S B S -> 4>X) and P{B S B S -> 4>£X) that 
a B S B S pair will yield a (f> or </>-lepton event. To esti- 
mate P(B S B S — > (j>X) we use the current world averages 
[1] of the inclusive branching fractions B(B S — > D S X), 
B(D S -> <j>X), and B{D -> <j>X). Here and in the fol- 
lowing D refers to the sum of and D° contributions. 
Also needed are estimates of the unmeasured branching 
fractions B(B S — > cc0) and B(B S -»■ DD S X). The former 
quantity accounts for direct B s —¥ <fi production, a sub- 
stantial fraction of which arises from B s to charmonium 
decays. We use the central value from the simulation, 
1.7%, which is roughly consistent with charmonium pro- 
duction in the B system. For the latter quantity we use 
a naive quark model prediction of 15% for b —> ccs. 

The inclusive <f> yield in B s decays can be expressed as: 

P(B S -> <(>X) = B(B S -> D&X) B(D S -> $X) 

+ B(B S -> cc</>) (6) 
+ B{B S -> DD S X) B{D ->• <j)X), 

from which we determine 

P{B S B S -> 0X) = 2P(B S -> - P(S S -> (7) 



The unknown quantities in Eqs. (2) and (3) are / s and 
the common normalization Rb- The ratio f s can be de- 
termined as a function of -Ecm by eliminating Rb be- 
tween the two equations. The result is presented in Fig. 
3. The ratio f a peaks around the T(55) mass. The to- 
tal excess below the B S B S threshold and deficit above 11 
GeV are consistent with zero within 1.5 and 1.3 standard 
deviations, respectively. 

Using Eq. (4), a \ 2 is constructed from the measured 
and expected values of P(B S B S — >• <j>£X) across the entire 
scan. The \ 2 is minimized with respect to B(B S £vX). 
The following processes contribute to C^e from B S B S 
events: primary leptons originating from a B s scmilcp- 
tonic decay, secondary leptons resulting from semilep- 
tonic decays of charmed mesons, and ^ or K ± misiden- 
tified as e ± or fi^. The contribution from primary lep- 
tons arises from events where one or both B s mesons 
decay semileptonically, and we determine the </>-lepton 
efficiency for each case (denoted for one semileptonic 
decay and e^ ee for two). It is found that ranges from 
8.5% - 10% and is about 10%. 




40 50 
Br(B -» I X) (%) 



FIG. 4. x 2 formed from the measured and expected yields, 
as described in the text, as a function of the semileptonic 
branching fraction. Note that since we express the branching 
fraction as the average of the e and y, channels, the physical 
bound is 50%. 

For the secondary lepton contribution, we consider 
events with up to two leptons coming from , D° or 
Df- decays. The selection efficiency in this case is esti- 
mated as the product of the <fi reconstruction efficiency in 
B S B S events in which neither B s decays semileptonically 
but a lepton candidate is identified (referred to below as 

), and a lepton detection efficiency determined from 
simulation (ef). It is found that lies in the range 
15% - 16.5%, and ef in the range 8% - 9.5% per lepton. 
The contribution from hadrons that are misidentified as 
leptons is estimated from simulation to be 3.3% of the 
(j)- lepton candidates in B S B S events. 

For the expected and measured 4> yields, we find: 
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e% t P{B s B s -»• <j>lX)? T -„ = (2e^ - e% u )B{D s 

+ B(B S £i/X)e5 



[-2 + S(P> S -> [£(P S -> &/X)] 2 
B(D S -»• </>X) + [1 - £(£ s -> ^Q]P(P S </>X) 



(8) 



e %lP{BsB s -> (^X) Secondary — 2e 



B(£) a -> + S(£> s -> £isX)B(D s -> <f>X) 



- B{D S -> tv<f>)B(D a -> [B(P S -> £z/X)] 2 

+ |p(P s -> <j)X){B{D s -> - 23(1), -> £z/X)) - £(P S -> 

- -> D S X)B{D S -> - £(P S -> D S X)B{D S -> lvX)B(D a -> 0X) 
+ B(P S -> D S X)B{D S -> tv<j))B{D s -> 

+ (B(D a -»■ </>X) - 2) ^ B(B,->DW-D j (X))B(D 

+ i3(B s -> D S X)P(B S [B(£>„ -> &/X) - £(£>« - 

+ i3(P s -> D S X)B(D S -> ^0) 

+ [£(P S -> + £(£> s </>X) - P(P S -»• (j)X)B{D s - 



-> &/-X") £(P S 



(9) 



£%i.P{B s B 8 
ei e P(B s B s 



■ (j>£X) Expecte d 
<f>£X) M eaauied = (1 - 0.033) ( C, 



x ^ B(B s ^dW-A(X))B(A^M) , 

i£u,d.s ) 

{e% e x 0.591 x S(P S ->■ - (2e^ - e^ £ ) x 0.289 x (fi(P 5 -> 

+e^ef [0.1375 - 0.2721 x £(P S -> &/X) + 0.1339 x (B(P S -> ^X)) 2 ] } 
/ s e£ + (1 - f s )e h (1 - f s )e <pl P{BB 



fsCh 



fs 



(10) 

(11) 



where that Eq. (10) is the sum of Eqs. (8) and (9) after 
substituting the values of known quantities. The first line 
in Eq. (10) expresses the contribution from primary lep- 
tons and the second that from secondary leptons. Eqs. 
(10) and (11) describe the measured and expected values 
used to form the \ 2 to be minimized, along with the sta- 
tistical uncertainties of each of the measured quantities 
and the uncertainties in the energy-dependent efficien- 
cies. 

The expression in Eq. (10) for the expected value of 
e^gP(B s B s — > 4>£X) is quadratic in the unknown B(B S — > 
£i>X), and so a x 2 formed from the deviation of the ex- 
pected from the measured values, summed over all bins 
above the B S B S threshold, is quartic in this unknown. 
Minimizing the x 2 with respect to the B s semileptonic 
branching fraction we find B{B S ->■ IvX) = 9.5+ 2 ;[j%. 
Figure 4 shows the dependence of % 2 on B(B S — > ivX). 

Systematic uncertainties are summarized in Table I 
and include the contributions described below. 

• Uncertainties for branching fractions, which are ei- 
ther taken from Rcf. [1] when known, or assumed 
to be 50% for B(B S -> cc(f>) and B(B S -> DD S X). 
These are separately listed in Table I, as is B(B S — ► 
D a X), which contributes a very large uncertainty 
compared to the other branching fractions. 

• Requirements used in the event preselection, in- 



cluding the lepton momentum requirement. The 
uncertainty due to the lepton momentum require- 
ment dominates in this group, and reflects the de- 
pendence of the result on the decay model used to 
simulate B a semileptonic decays. 
Fixed parameters used in the fits to tukKi includ- 
ing to^, T^, cr. 

The parameterization of the background and ab- 
sence of a term in the fit corresponding to threshold 
contributions from light scalars. 
Uncertainties in particle identification (PID) effi- 
ciencies and hadron misidcntification probabilities. 
The determination of P(BB -> <j)X) and P(PP -> 
<j)lX) in Y(4S) data (these quantities are deter- 
mined to 1% and 1.8% relative uncertainty). 
Sensitivity of efficiencies to differences in branching 
fractions implemented in simulation compared to 
their measured values. 

Uncertainties in the continuum-subtracted number 
of events due to ISR and two photon events, which 
do not follow a 1/Eq M energy dependence. 
A correction made to the continuum subtraction 
of the number of PS-like events due to an over- 
subtraction found in simulation studies. The size 
of this correction is about 1% of the amount to be 
subtracted; we use ±100% of this correction as a 
systematic uncertainty. 
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TABLE I. Relative multiplicative and additive systematic un- 
certainties for the measurement of B(B 3 — > £vX). 



Multiplicative Systematics Relative 

Uncertainty (%) 



B{B S ->■ D^X) 


+8.72/- 13.58 


B(B S — > cc</>) (Unmeasured) 


±3.20 


S(B S ->■ DD S X) (Unmeasured) 


+1.12/— 1.16 


Other Branching Fractions 


+0.52/-0.54 


Event and Lepton Selection 


+1.99/-2.85 


Fixed Fit Parameters 


+0.49/-0.15 


Background Parameterization 


±0.93 


P1D and Lepton take Rate 


±3.21 


P(B u , d B Utd -»■ 0) 


±1.47/— 1.69 


Simulation Branching Fractions 


±2.09 


ISR and 27 Background 


+1.57/-7.14 


Correction to Event Subtraction 


+1.88/-4.59 


Technique bias 


+0.39/- 10.00 


Total Multiplicative 


(±10.87/-19.92)% 


Additive Systematics 


Uncertainty (xl0~ 3 ) 


Other Branching Fractions 


+0.56/-0.64 


P(B u , d B u , d -> 4>lv) 


+4.30/-3.90 


Total Additive 


(+4.34/-3.95) x 10~ 3 



Total Systematic (±11.20/-19.34)xl0 -3 



In conclusion, we performed a simultaneous measure- 
ment of the B s semileptonic branching fraction and its 
production rates in the CM energy region from 10.56 GeV 
to 11.20 GeV. The semileptonic branching fraction is con- 
sistent with theoretical calculations in Refs. [9] and [10]. 
Our measurement of the B s production rates are con- 
sistent with the predictions of coupled channel models 
[11], in which B s production peaks near the T(5S) and 
is vanishingly small elsewhere. 
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• Possible bias in the x 2 minimization technique at 
low statistics. Firstly, evaluating the behavior of 
this method for extracting B(B S — > IvX) for many 
pseudo-data samples derived from the simulated 
dataset gives evidence for a small bias at low statis- 
tics. Secondly, it was found that the analysis per- 
formed in high statistics simulation tends to overes- 
timate B(B S — >• ivX) by an amount corresponding 
to half the statistical error reported. 

To determine whether the uncertainties from these 
sources scale with the result or not, each was evaluated in 
a simulation sample with a higher semileptonic branching 
fraction and compared with the result in the normal sim- 
ulation sample. It was found that the uncertainty from 
the determination of P(BB — > <p£X) in T(4S I ) data does 
not scale with the branching fraction, nor does the uncer- 
tainty contributed by several of the input branching frac- 
tions. These are thus separated in Table I. The remain- 
ing uncertainties are found to scale with B(B S — > IvX) 
and thus to be multiplicative. 

Our final result for the inclusive semileptonic branch- 
ing fraction is 9.5l" 2 /o-i.9%' which is the average of the 
branching fractions to e and fi. 
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